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Mean Field Theory has been extensively used in the study of systems of anyons 
in two spatial dimensions. In this paper we study the physical grounds for the 
validity of this approximation by considering the Quantum Mechanical scattering of 
a charged particle from a two dimensional array of magnetic flux tubes. The flux 
tubes are arranged on a regular lattice which is infinitely long in the "y" direction 
but which has a (small) finite number of columns in the "x" direction. Their physical 
size is assumed to be inhnitesimally small. We develop a method for computing the 
scattering angle as well as the reflection and transmission coefficients to lowest order 
in the Aharonov-Bohm interaction. The results of our calculation are compared to 
the scattering of the same particle from a region of constant magnetic field whose 
magnitude is equal to the mean field of all the flux tubes. For an incident plane wave, 
the Mean Field approximation is shown to be valid provided the flux in each tube is 
much less than a single flux quantum. This is precisely the regime in which Mean 
Field Theory for anyons is expected to be valid. When the flux per tube becomes of 
order 1, Mean Field Theory is no longer valid. 
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I. INTRODUCTION 



There has been much interest in recent years in the phenomenon of fractional statistics in 
two spatial dimensions. Particles possessing fractional statistics, known as anyons, besides 
being interesting in their own right, have found applications in the Quantum Hall Effect 
and have also been investigated in the context of High Tc Superconductivity. 

Anyons are often studied in the framework of "Chern-Simons" theory in which they 
are modeled by "attaching" an infinitesimal tube of "statistical" flux to each particle in 
the system. In order to calculate various properties of the anyon gas it is often assumed 
that each anyon "feels" as if it were traveling in the mean magnetic field due to all of 
the other anyons in the system. This "Mean Field Theory" has been quite successful at 
describing many properties of anyon systems including Anyon Superconductivity and 
the Fractional Quantum Hall Effect 0. In fact in Chern-Simons Field Theory, the mean 
field approximation can be shown to be valid when the statistical flux per particle is a 
small fraction of a flux quantum. Classically this Mean Field (MF) idea seems absurd 
since the magnetic field is zero except at isolated singularities and so there is no Lorentz 
force. Quantum mechanically however the vector potentials themselves attain importance 
@]. Thus it seems at least plausible that a particle moving in a region populated by flux 
tubes would in some respects behave as though it were in a mean magnetic field. It is this 
idea which we investigate in this paper. 

It is quite straightforward to see that Mean Field Theory is at best valid only when the 
flux per particle is not too large. Assume that a given anyon "feels" a mean field and hence 
travels in a circular (Landau) orbit. If the area enclosed by the orbit contains, on average, 
many other anyons then the approximation may be valid. Proceeding in this way one finds 
both for models based on bosons |^ and on fermions that the MF approximation is 
self-consistent only when the flux per particle is much less than a single flux quantum. 

It is our goal in this paper to examine the physical ideas implicit in the above argument 
in more detail. One question of interest in this context is: To what extent does a particle 
traveling in the presence of a number of infinitesimally small flux tubes really behave as 
though it were in some mean magnetic fleld? In other words, to what extent does the 
presence of the gauge fleld mimic the effect of a spatially constant statistical magnetic fleld? 
Ideally we could consider an array of inflnitesimally small flux tubes and study the quantum 
mechanical motion of either a charged particle or another flux tube in the presence of this 
array. This however turns out to be a very difficult problem to solve 0. We choose instead 
to consider a specific configuration of fiux tubes consisting of a lattice which is of infinite 
extent in the y direction but which contains a (small) finite number of columns in the x 
direction. A related problem (in which the array is random) has been considered previously 
in Ref. in which scattering from the array was considered as an incoherent sum of the 
scattering off each flux tube. In this paper we consider the coherent scattering of a charged 
particle incident on a lattice of flux tubes. The method which we develop is a variation on 
flrst order scattering theory. Our main objective will be to compute the scattering angle 
and to compare the result with scattering from a strip of constant magnetic fleld. 

We thus begin in Section II by calculating the scattering angle of a plane wave incident 
on an inflnitely long strip of constant magnetic fleld. In Section HI , after calculating some 
analytic results, we develop a formalism, similar in spirit to the approach conventionally 
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used to compute the index of refraction from the microscopic scattering amphtude P,p!0| , to 
study scattering of a charged particle from a lattice of flux tubes. We then use this formalism 
to find an expression for the scattering angle to first order in the interaction, which we then 
compare to the analytic results. (Note that while exact solutions are available for scattering 
from one ^ and, more recently, two ||rT| flux tubes, the problem of scattering from an 
infinite array of tubes cannot be completely solved analytically. We will find however that it 
is possible to solve for the scattering angle exactly in a limited number of cases.) In Sections 
II and III we also discuss the reflection and transmission amplitudes. Section IV contains a 
discussion of the results. 
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II. SCATTERING FROM A STRIP OF CONSTANT MAGNETIC FIELD 

In this section we briefly review the result of scattering a charged particle from an 
infinitely long magnetic strip in two spatial dimensions. We thus consider a situation in 
which the magnetic field B is constant in the region < x < d. It is convenient to work in 
a gauge in which Ar^ = and 

' -Bd/2, if X < 

Bx-Bd/2, ifO<x<d (1) 

Bd/2, ifx>d 

The Schrodinger equation for a particle with charge e subject to the above magnetic field is 
given by 

- ^[V - ieA]^^P{x,y) = Eij{x,y). (2) 

This equation is solved by separating the x and y dependence and representing the latter as 
a plane wave. We thus write 

ij{x,y) =e'^''y+'^^yf{x) (3) 

with 

-1 j2 jdj 

Suppose that the scattering particle is incident from the right (positive x). The incident 

wave 

i^inc = e*('=«+^)3/e-''=-" (5) 

will solve the Schrodinger equation for x > d provided 2mE = k1 + k"^. 

In the interior region (0 < x < c/) the particle behaves as if it were in a potential which 
is a truncated parabola. Although an explicit solution is in general not available, it is rather 
straightforward to obtain a series solution. For certain combinations of the parameters the 
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series truncates and the wave function is simply a Hermite polynomial but in general the 
series is infinite. (In the case when the series truncates, the wave function is related to that 
of the Landau level.) For x<0 the wave function is again given by a plane wave 



fix) 



but now 



It thus follows that for x<0 



2mE = kl + {ky + eEdf 



= Jkl - 62^2^2 _ 2kyeBd 



(6) 



(7) 



(8) 



Transmission of the particle to the region x<Q can occur if the energy is sufficiently large 
so that A;2 > e^B'^d? + 2kyeBd. In this case the solution to Eq. (0) can be written as 



r e 



-^V2 (coGo(X) + ciGi(X)) 



X <0 
0<x <d 
X > d 



(9) 



where k^^ is given by Eq. and 

X = J\eB\ I X 



(10) 



is proportional to the distance from the centre of the parabola. The functions Gq and Gi 
are given by 



G'i(X) =X + 



2(1-A)^3 , 22(3-A)(l-A) 



3! 



5! 



(11) 
(12) 



w.th A4(^|tS_,J (13) 

Co, ci, p, and r are constants which are fixed by requiring that the wave function and 
its derivative be continuous across the two boundaries, p and r are the refiection and 
transmission amplitudes respectively. Note that the series ( pH]) and (0) converge for all 
X. (To see this note that for the high order terms in the expansion the ratio of successive 
terms behaves like the corresponding ratio in the expansion of exp(2X2) |T^.) If the energy 
is small so that k^. < e^B'^d^ + 2kyeBd there is no transmission and one obtains instead a 
decaying exponential for x < so that the wave is totally refiected. 
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It is quite straightforward to compute the angle by which the particle is scattered. The 
simplest way to do this is to compute the gauge invariant current density given by 

J = — Im [i/j* (V - ieA) tp] (14) 
m 

The angle (j) which the current makes with the x-axis is then given by 

tan0=-^. (15) 

Using Eqs. (|]) and (||) we can compute the angle of both the incident and the transmitted 
wave. We find that 

k 

tan (0inc) = 7^ 

tan (0trans) = , = (16) 

Ikl - e^B^d? - 2kyeBd 

This agrees exactly with the result obtained for the scattering angle in a purely classical 
treatment. In fact the result simplifies significantly when \eB\d << k^ in which case 

, , , , eBd . 

tan (0trans " 0inc) - -7— (17) 

From Eq. ([T6D we can understand the physical origin of the restriction which was imposed 
above for k^. When k^ = e^B'^d^ + 2kyeBd the angle of the transmitted wave 0trans = ±7r/2 
so that for k^ < e^B'^d'^ + 2kyeBd the particle "turns around" and is re-emitted on the right. 

It is evident from the above discussion that in order to find 0trans given kx and ky it is 
not necessary to know the detailed form of the wave function in the interaction region. In 
fact all we actually used to get the scattering angle in the above analysis was our knowledge 
of the vector potential on both sides of the barrier and the fact that ky was conserved. We 
shall come back to this point later. 

We now proceed to a computation of the reflection and transmission amplitudes p and 
r. In general exact values can only be obtained numerically since the expressions will 
necessarily contain the (generally infinite) series G'o(X) and Gi(X) given in Eqs. (0) and 
(|l^) and their derivatives. It is however possible to obtain approximate analytic expressions 
for p and r when the magnetic field is small by first relating G'o(X) and Gi(X) to the 
confluent hypergeometric functions as follows: 

G'o(X) = M(-^,i X^), (18) 

G'i(X)=XM(l-^,|x2), (19) 



where 



M a, 6, ^ = 1 + + + . . . + + . . . , 20 

(0)22! (O)n^! 

(a)„ = a(a + l)(a + 2)---(a + n-l), (21) 

(a)o = 1. (22) 
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The confluent hypergeometric functions are in turn related to the Bessel functions by 



1 oo 

M(a, b, z) = r(6)e^(-te - az)"^-"^' ^ ^n(-)^(& - 2a)-t J,_i+„((2^6 - 4za)^), (23) 

n=0 ^ 



where -Ro=l) -Ri=0, R2=b/2, and 

i?„+i = [(n + 6 - l)i?„_i + (2a - 6)i?„-2]/(ri + 1), (24) 



for n>2 13 



We shall compute p and r for small eSc?^ only in the case ky=0. To do this we shall 
need the values of G'o(X), G'i(X) and their derivatives at both x=0 and x=d. When x=d, 
X=0, so these values are easily found by direct substitution into ( pTf ) and ([T2|) and the 
derivatives of these sums. When x=0, however, X=—^\eB\d, so we use (p!8|), (p!9D, and ( |23D 

to obtain expansions for Gq, Gi and their derivatives in powers of the quantity ^\eB\d. The 
calculations are considerably simplifled by the fact that "6" in ( [23| ) is always a half integer so 
that the various Bessel functions reduce to simple trigonometric expressions. Furthermore 
the arguments of the Bessel functions are simply equal to "fcrf". After a rather lengthy 
calculation we find, for \eB\d'^«l and k»\eB\d 

e^B'^d'^ ( siv?kd .,sm2kd 1 A ^, 0^0,= , 

Note that the above result satisfies the conservation of flux 

\p? + t\^? = 1 (27) 



V 



to the required order, 0[e B d ). 



III. SCATTERING FROM AN INFINITE STRIP OF FLUX TUBES 

A. Analytic Results 

In this section we consider the scattering of a charged particle from a regular array 
of inflnitesimally small flux tubes. Our goal is to establish conditions under which this 
scattering is equivalent to the scattering from a uniform magnetic fleld whose magnitude is 
the mean fleld of all the flux tubes in the array. 

To this end consider a situation in which we have an array of inflnitesimally small flux 
tubes in two spatial dimensions each with a flux $. We take the lattice to be inflnite in 
extent in the y direction with the flux tubes separated by a distance ^, but to have + 1 
"columns" in the x direction with flux tubes positioned at x = 0, ^, . . . , NC,. We shall 
often choose to have only one column of flux tubes (A^=0). (The fact that we have set the 
spatial separation of the flux tubes in the x and in the y directions to be equal is simply a 
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convenience. Our results are easily generalizable to the case of unequal spacing.) We choose 
the gauge potential of each individual flux tube in a radial gauge relative to the location of 
the flux tube so that 

Ar = (28) 



Ae = — = -^ (29) 
iTcr er 

where r and 6 are the coordinates relative to the position of the flux tube and ( provides a 
convenient parameterization of the flux. 

The contributions from each of the flux tubes in the lattice can be summed up to obtain 
a closed analytic result for the full gauge potential. To do this one begins with the result 



1 °° 1 

cotvrz = — , (30) 

for complex z. From this one can derive the formulae 

^ [ma — pY + Y 2a;7 a a 



m=— oo 



and 



E , "^"T/ =-^[cotvr(^)+cotvr(^)] (32) 
^ (ma-/5)2 + 72 2a^ ^ a ' ^ a ^ ' 



m=~oo 



where a, (3 and 7 are real. Applying these results to our problem gives 

(vr ^ sin(^) 
^^ = Vi to cosh(f(x-<))-cos(^) ^^^^ 

A - Cvrf sinh(f(x-<)) 

ee„^ocosh(f(x-0)-cos(^) ^ > 

where the sum is over the iV + 1 columns of the flux tube array. In fact when A^=0 the sum 
collapses to a single term. 

Asymptotically (for large positive and negative x) the gauge potential becomes 



^ as X — >■ ±00 

(35) 



Ay^^^{N+l) as x^±oo 



This asymptotic value is actually reached quite quickly. In fact the first order corrections 
damp exponentially as exp(-27r|x|/^). Not surprisingly the asymptotic value of the field is 
precisely that of the constant magnetic field B given in the previous section (Eq. (||)) with 
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B - ^ (36) 

and 

d={N + l)i (37) 

so that B can be interpreted as the mean field of the fiux tubes. 

One of our goals in this section is to compute the angle by which a charged particle 
incident from x = +00 is scattered by the array of fiux tubes and to compare the result 
with the scattering angle off a strip of constant magnetic field. We thus wish to compute 
the transmitted angle 0trans given kx and ky which can be defined asymptotically as in 
the previous section. In the case of a strip of constant magnetic field we found (in the 
previous section) that in order to compute the scattering angle it was possible to ignore 
the behaviour of the vector potential and the wave function in the region 0<x<d and to 
simply use the asymptotic forms of these quantities and the fact that ky was conserved. It 
is tempting to suggest that the same result should apply in the present situation, in which 
case the scattering angle for the array of fiux tubes would be the same as that for the strip 
of constant field. This result is however only true under certain restricted conditions. The 
technical problem with the argument is that in our case, due to the lack of (continuous) 
translational invariance in the y direction, the Schrodinger equation does not separate into 
X and y pieces. The value of ky is thus not conserved through the "interaction" region of the 
potential. The discrete translational invariance which is present in our case does, however, 
lead to interesting predictions. This will be discussed below. 

The simplest way to see that the above argument would lead to false conclusions is to 
consider the case when the fiux in each fiux tube is an integer. In this case the corrections to 
the "Mean Field" result are large since the beam is in fact totally undefiected. This follows 
from the fact that the interaction can be completely "gauged away" by a nonsingular (except 
possibly at the centres of the fiux tubes) single-valued gauge transformation. In fact when 
( is an integer the wave function 



^ = e^^-'^-^+'^^^^exp 



i( ^ tan 



iCiry 



(38) 



n=0 \ — 

is single-valued and solves Eq. (^ with A given by Eqs. (|33|) and (|3^) . The resulting 
current has tan0 = ky/kx everywhere (except possibly at the singularities). In cases where 
C, is not an integer the wave function (|38D is multi-valued and is therefore not an admissible 
solution to the Schrodinger equation. These results were, of course, already noted in the 
original paper of Aharonov and Bohm . It is thus clear that replacing the fiux tubes by a 
mean field will in general give an incorrect answer. We shall study below the circumstances 
under which the mean field result is correct. 

As one might expect from the fact that there is no scattering when C, is an integer, the 
current (and hence the defiection angle) is a periodic function of C- To see this, suppose i/'c 
is an exact solution to the Schrodinger equation with < C < 1. Let p be an integer and 
perform a gauge transformation with gauge function 

2iTy\\ 

X = - 2. tan ( —72^, ^ I - (39) 
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This now gives an admissible solution for the case ( P + C- The current is of course 
unchanged by the gauge transformation. Thus the scattering angle is periodic in ( with a 
unit period. 

It turns out that for certain values of kx, ky, and ( we can solve for the scattering angle 
exactly without knowing the form of the wavefunction in the scattering region [|T5|. The 



solution depends on the periodicity of the lattice in the y direction |]I6[. Following the 
usual discussion of Bloch's theorem |T^, one can define a unitary translation operator as 
follows: 

TJ{x,y) = f{x,y + nO. (40) 

Since our Hamiltonian (with the vector potential defined in (|33| ) and ( p^ ) is invariant under 
y ^ y + n^, we have the operator identity 

[Tn,H] = 0. (41) 

Furthermore, the T„'s form a representation of the group of translations by so that 

T-nTfyi TjYiT^ Tn+m,} (42) 

which means that eigenstates of H may be chosen to be simultaneous eigenstates of all of 
the Tn- 

Hifj = E-^, (43) 
Tni^ = Cni^, (44) 

with CnCm=Cn+m- Thus c„=exp(27rm(5) for some real parameter S, which will be related to 
ky in this case. 

We now consider the form of the wavefunction as x —>■ oo (the incoming wave). In this 
region we may take (see Eq. (|]) ) 

tP{x, y) = e^('=^-C-(^+i)/?)f (e-^^-^ + pe''^-"), a; ^ oo, (45) 

so that 

S.td.<A!l±}l, (46) 

To the far left of the scattering region we again seek a plane wave solution, so we set 

^(x, y) = re''^k'y-'i'''^^+^)/Oy^-i~k.x^ ^ _^ (47) 

The fact that c„=exp(27rm(5) then gives 

,/ , 27rm 

K = ky + ^- (48) 
Energy conservation then determines k^'- 
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97r dirk 1 

= [kl - {jfiCiN + 1) - mf + ^(C(iV + 1) - m)]^ (49) 

so that if is real the angle of the beam of transmitted particles is determined by the 
equation 

ky - ^{C{N + I) - m) 
tanCans = -0 ^^^^^ TT r • (50) 



[kl - (WiCiN + 1) - my + ^(C(iV + 1) - mp 



Comparing the above equation with Eqs. ([T6|), ( p6D and (|37D, we see that agreement with 
the MF result is achieved when m=0. 

We now restrict ourselves to the case where N=0 and consider plotting tan</)trans versus 
(. The "critical points" (i.e., the points for which k^=0 ) are given by 



C;^ = -^ + |7±(^^ + ^J)^|^- (51) 



If {kl + ky)^ ^<Tr, the scattering angle is well defined in the intervals {(~, Cm) is given by 
Eq. (|50D . In the intervals(Cm! Cm+i)) the k^ are imaginary and the beam of particles is 
backscattered by the flux tubes. If {k^ + fc^)2^>7r, there are regions of the ( axis for which 
more than one scattering angle could be defined. Had the incident beam been characterized 
by a wave packet in our calculations, the scattered particles would have split (asymptotically) 
into several wave packets with distinct directions in space, given by (|5(]|). Since our incident 
beam is characterized by a plane wave the scattered wave can never be separated into its 
various pieces and the scattering angle is not well defined. If we view the situation from the 
perspective of diffraction theory, we see that the "circular waves" being emitted from each 
of the scattering centers are interfering constructively in more than one direction. 

A similar result holds for any value of N. There may be certain values of the parameters 
for which there is no value of m with real /c™. In such cases the beam of particles is totally 
reflected by the flux tubes. If there is only one integer m for which k^ is real the scattering 
angle is well defined and is given by (|50|) . If, in addition, m=0, the angle agrees exactly 
with the MF case. Finally, it is possible that there exist more than one integer m for which 
fc™ is real. This situation corresponds to the "diffraction" case which was considered above. 

Finally, note that the invariance of 0trans under ( ^ ( + m (which is expected in light 
of the gauge invariance arguments given above) is already included in Eq. (|50|) by simply 
taking m — >• m — m (iV + 1). 



B. Numerical Results 

The above analytic analysis is somewhat limited in that it only predicts the scattering 
angles and in that it depends in an essential way on the periodicity of the lattice. We shall 
now develop a numerical technique which will allow us to examine not only the scattering 
angle but also the lowest order values of the reflection and transmission amplitudes. This 
technique can be generalized to a non-periodic array of flux tubes. 

As discussed in the introduction, we expect the Mean Field result to be valid when the 
flux per "particle" ( is small. We thus analyze our problem for small (. Consider a wave 
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which is incident from x = +00. This incident wave will induce a scattered wave from each 
of the flux tubes in the lattice. To linear order in ( the total wave function can be viewed as 
the incident wave plus the sum of the scattered waves from each individual flux tube. This 
procedure is, however, significantly more complicated than in ordinary scattering theory due 
to the presence of the Aharonov-Bohm phases. 

In order to proceed we shall begin by reviewing the results for scattering from a single 
flux tube which we assume is positioned at the origin. In this case, an exact scattering 
solution is known 0]. It is given by 

00 

i;= Y: e-^l™+^lJ|„+c|(^Oe^"^^ (52) 

m=— 00 

where J is a Bessel function and the angle 6 should be taken from — tt to vr (where ^=0 
corresponds to the positive x axis). To see that this is a scattering solution one notes that 
as X — s> 00 

so that the incident current density is in the negative x direction. In fact for \(\ < 1, 
Aharonov and Bohm (with modiflcations due to Hagen \TE\) showed that ( |52D reduces asymp- 
totically (as r 00) to 

pikr 

*^^^""""""^'(2SS)^^'"^'^SW2)' '''' 

where the upper (lower) sign holds for Q positive (negative). The above asymptotic expres- 
sion for ip requires one further crucial assumption which is that {kr{l ± cos6'))^/^ are large. 
This assumption fails in the forward and backward directions so that ( ^41) is inaccurate in 
these regions. This is fortunate since (^) is convergent and single-valued for all 9 as can be 
conflrmed by a direct evaluation of ( |52D near ^=±7r. The asymptotic form (Q), however, is 
not single valued and it diverges when ^=±7r. 

This "splitting" of the wave, in the case of a single scatterer, into an incident and a 
scattered wave is essential for us since it is the starting point for summing up the scattered 
waves from a lattice of flux tubes. We cannot simply use Eq. ( |5^ for this splitting, even 
though the "scattered" wave is proportional to sin((^7r) which is small for small since Eq. 
(p^) is not valid for all 6. Eq. (|^) does however suggest a more appropriate way to deflne 
the scattered wave. We consider the deflnition: 

^scat = ^fuu - e-^'^^-^^ (55) 

where ipf-^w is the full wave function given in Eq. (^2|) which can be evaluated (for example) 
numerically. This deflnition is not without problems since, of course, ipscat is not single- 
valued. There are however several ways to circumvent this problem and to obtain useful 
results with this approach. 

The most aesthetically pleasing approach to dealing with this "multi-valuedness" prob- 
lem is to work with the currents rather than with the wave functions whenever possible. We 
begin by taking the covariant derivative of (^21) directly. This gives the simple result 
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(V - ieA)ij = -nkij + if, + jfy, (56) 

where 

= |[e^^^/'(e-^V_c+i - ^J-c) + e-^<-'\e~^' J^^, + zJ^)], (57) 

/. = ^ie'^^-'^ie-^'J-^+i - J-c) + e-«-/2(ze-^ Jc-i + ^c)]- (58) 

The Schrodinger equation can now be used to derive an equation satisfied by fx and fy. 
Applying the covariant derivative to ( p6D gives the relation 



df df 

- ^^/- + ^ + ^ - ^e^-/- - ^^^yfy = 0' (59) 

which may be verified explicitly. 

For the case of a lattice of fiux tubes the covariant derivative of the total wave function 
may be obtained, to leading order in (, by adding, to the term —iikiptotui, the sum of the 
contributions of the f^^s and /^'s from each fiux tube. We shall specialize to the case when 
there is only one column of fiux tubes and when ky=0 (i.e. normal incidence). We shall 
discuss the more general case later. Let us define the vector n{x,y) as the sum of the 
contributions of fx and fy from all the fiux tubes in the lattice. That is 

oo 

^x{x,y) = fxikrn,9n,C), 
n=—OD 
oo 

ny{x,y)= Y fyikr^,Or„C), (60) 

n=— oo 

where 

rn = {x' + {y-nO'y/', (61) 

and where 6'„ = vr is parallel to the x-axis and corresponds to the direction of propagation 
of the incident wave. To lowest order in ( the total wave function satisfies 

(V-ieA)^ = (-iiA;)^ + ri (62) 

This is analagous to what is done in ordinary scattering theory where the scattered waves 
from each individual atom are summed to obtain the total scattered wave. 

To further check the self consistency of the above procedure we apply the covariant 
derivative to Eq. (^). \E' is expected to satisfy the Schrodinger equation far away from any 
of the fiux tubes. (In fact in our case it will turn our that "far away" may be as close as a 
few lattice spacings.) It will do so provided O satisfies the equation 

dQx dQy 

— ikVt^ H — 1 — — ieAx^x — ieAyVly = (63) 

dx dy ^ ^ V ^ 



where Ax and Ay are the total gauge potentials given in ( |55| ) and (0). Consider the 
asymptotic region on the left where Eq. (|35|) holds. Applying (|59D to every term in the 
sums in (Bl) gives 
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E 



Tcrt 



rn 

' J X 




0{C 



TT 



(64) 



It is easy to see that the /'s are, asymptotically, of order (it. Thus the Schrodinger equation 
will be satisfied asymptotically provided fl is also of order ^vr. Eq. (^) guarantees, that 
this will be the case. 

We have not been able to obtain an analytic expression for ft, even asymptotically. 
We have thus chosen to evaluate O numerically. Some details of the numerical work are 
discussed in the Appendix. We have found that to a very high precision and for a range of 
k which is also discussed in the Appendix (basically 2t[ > kS, » 2^7?) 



2i 



n ~ --sin(C7r) e-^'=^j 



(65) 



This result is valid provided one is not too close to the array of fiux tubes. Note also that 
to this order in C, the value of k (which is actually kx) is the same for the incident and the 
transmitted wave. Also, as expected, when C^n is small (|^) is satisfied. 

The next step in evaluating the scattering angle is to compute the current. Using (|62[) 
the current is found to be 



J = — Im[**(V-zeA)^], 
m 

1^ .-.^,2 



+ Im(^*ri)] 



(66) 
(67) 



This expression could potentially cause us some difficulty since it involves ^ directly whereas 
we have only computed the covariant derivative of Fortunately, to the order in to which 
we are working, there is no ambiguity in \1/ and we can take \l/ to be simply the incident 
wave since linear order corrections will lead to quadratic order corrections to the scattering 
angle. Thus 



-ikx 



+ 0{Cn) 



and 



= l + 0(Cvr) 



We are now ready to use Eq. ( |67D to evaluate the scattering angle 0trans- 



tan0 



trans 



+ 0(CV) 



(68) 



(69) 



(70) 



Eqs. (|T^) and (|36D can now be used to compare this result with the scattering off of a strip 
of constant magnetic field. We see that to lowest order in (tt the scattering from an array of 
fiux tubes agrees exactly with the Mean Field result (which coincides with the analytic result 
Eq. (|50D in this case). This agreement provides a good check of our numerical method. 

In the previous discussion we have avoided summing the non single-valued scattered 
wave functions by summing only their covariant derivatives. It is however possible to sum 
the scattered waves directly without first evaluating Dip. Using the definition (|55|) of the 
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scattered wave from an individual flux tube we note that the scattered wave at {x, y) due 
to a flux tube located at is given by 

where ip is the full wave function for each individual scatterer and r„ and 6n are defined as 
in the previous discussion. To lowest order in ( the total scattered wave will be the sum of 
the waves scattered from each individual flux tube: 

oo 

*scat= E Ccat(fcr„,e„,C). (72) 
n=— oo 

We know of no way to do this sum analytically. We have thus evaluated (|7^) numerically 
by summing the terms until satisfactory convergence was achieved. For fcr„ large and On 
near ±7r/2 the asymptotic form ( |5^ ) for ip for each flux tube was used. In other cases the 
sum (^) was evaluated directly. This sum converges fairly rapidly for small and moderate 
values of A;r„ since for large index 

Mx) - (73) 

The numerical work showed that to the left of the scattering region and for momenta in 
the range 27i>kC,» 2(7t the scattered wave has the approximate form 

^scat ^ e-''^"(e(x, y) + zsm{Cn)x{y)). (74) 

For not too small e was found to be negligible compared to (it and x "was found to have 
the approximate form: 

X{y) ^ -^(^^ + for n^<y<{n + 1)^ (75) 

for oo < n < +00. In fact x does have some oscillatory x dependence but it was found to 
be negligible compared to the y dependence of x- The slope "—2/^" is accurate to about 
.2%, and (fc^ — l)/k^ is accurate to about 1%. For a more thorough discussion of e and x 
and of the limits on the momentum we again refer the reader to the Appendix. 

It is clear from our method that \l/scat must have discontinuities along lines starting at 
the site of each flux tube and extending out to the left parallel to the x axis. This is because 
we have subtracted out the multi-valued incident wave exp{—ikx — iC,6) from each (single- 
valued) wave function. This brings up the difficult problem of what to choose as the incident 
wave in the multi-vortex case. Although it is true that to lowest order it is simply a plane 
wave we must be very careful since the scattered wave (as we have defined it) is not single 
valued. 

One important condition which \l/inc should satisfy is 

(V - ieA)^inc = -nk^ inc (76) 

The form \l/inc =exp(— ifcx + iC'^y/^) satisfies this condition asymptotically (when x<<0). 
However due to the form of x which has discontinuities resulting from its multi-valuedness 
it is important to choose \l/inc to have the same discontinuities. In fact 
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^inc = e-'""^ {l + iCiT^y-^^ for ni<y<{n + 1)^ (77) 
also satisfies Eq. (^) asymptotically but only to order C,n. Setting 

^ = ^inc + ^scat, (78) 

we can compute the covariant derivative of the total wave function \E'. To linear order in C,% 

(V - ieA)^> = -iik"^ - 3jCrce-''''= (79) 



This is precisely the result which we obtained using our previous method (see Eqs. (|62D 
and (|65D ). We thus see that with this method as well, we recover the Mean Field result 
when the flux per particle is small. Note that we have shown that is a plane wave with 
corrections linear in (tt. This was assumed in our previous discussion (Eq. (pBD). 

It likely that the linear correction to the plane wave above should be exponentiated. 
From this point of view the discontinuities are very useful since they allow the perturbative 
corrections to ip to remain small even for large y. 

It is by now clear that the two approaches discussed above are complementary. In the 
former approach it appeared as though we could avoid any problems of discontinuities but 
in fact the discontinuities were all "hidden" in The latter approach has the advantage of 
showing clearly what is going on by evaluating the wave function directly. From a calcula- 
tional point of view the former approach is advantageous (at least in lowest order) since the 
derivative of \I' is found directly without having to evaluate at several points and perform 
a linear regression to find the slope. 

It is a simple matter to extend the above results to the case where ky is nonzero. A 
subtlety which one encounters numerically is that the wave functions (or their covariant 
derivatives) due to each flux tube must be multiplied by appropriate phases in order that the 
plane waves incident at each site are in phase. Our two methods again give complementary 
results so we describe here only the "first" method in which the covariant derivative or ^ is 
taken first. The appropriate generalization of (p^) is 



(V - leA)^ = {-iik^ + ]iky)<i' + fl, (80) 
where f2 has been studied numerically and is approximately given by 

n ~ -|s^n(C7^)e^('=-^-'=^^) fl^i + jl (81) 

in the asymptotic region on the left. We refer the reader to the Appendix for more details 
of the numerical results. Generalizing from the result for ky = Q we take 

^ = gi(M-fc.^) + 0(^7^) (82) 

We then find 

tan(0trans - 0inc) = + ^'(C'vr'). (83) 
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This result again agrees with the Mean Field result (Eq. ([TTD) and with the exact result, 
since "m"=0. 

The extension of our results to the case where there are several columns is also a simple 
matter provided the perturbative analysis remains valid. Suppose we were to add a new 
column positioned at x=C,- In the local coordinates (x , y ) of that column an incident plane 
wave exp{ikyy — ik^x ) gives rise to 

n' = J'sin{CT^)e'^'''y-^-^'^ f + j V (84) 



Since y =y and x =x — ^ the wave incident on the column at x=^ must be multiplied by the 
phase exp{—ikxC,) to obtain the correct phase of the plane wave at the locations of the two 
columns. Thus 

Qtot = O + e-'''^^n' = 20 (85) 

so that 

tan(0trans - 0inc) = + '^(C^^T^) (86) 



for A^=l. In general, for + 1 columns. 



tan(0,,,.3 - 0inc) = -^^^^^^^ + OiCV). (87) 

The net effect is that our perturbative expansion is now in "(A^ + l)(^7r" rather than in "(^tt". 
Thus even for multiple columns the scattering angle agrees with the exact result and thus 
with that obtained from the Mean Field analysis provided the flux per particle is sufficiently 
small. Again this agreement provides a good check of our numerical method. 

So far we have concentrated on comparing the scattering angle for the lattice of flux 
tubes with that obtained in the Mean Field analysis. We now describe a comparison of 
another quantity of interest namely the transmission and reflection amplitudes r and p, for 
which there are no exact results available. We shall discuss only the case ky=0 for which we 
shall compute numerically the sum in Eq. (ff^). 

In the previous section we computed p and r for the Mean Field case. In terms of the 
"flux parameters" ( and ^ and in the limit when (it is small and kC,»\C\7i they are given 
by (see Eqs (H, (|3), (HD and (H)): 



Pmf = e % TT + 1{ - -tt-^t) + 0{Q TT ), (88j 



V {kO' '2{kO' {kO 



x..^l + CV^^---j+0(CV). (89) 

(The extra factor of e*'^^ has been added to compensate for the positioning of the mean field 
relative to the fiux tubes.) Note that to order (tc we have the simple result that p=0 and 
r=l. To this order the wave is entirely transmitted in agreement with the classical result. 
We now describe the calculation for p and r for the lattice of fiux tubes. 
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We begin by summing ([72|) numerically. When k^»2\(\'n' and x>0 we find 



^refi = -^smCne""^, (90) 

to an accuracy of about 1% in the coefficient of "— i/Zc^". A naive interpretation of this 
result would imply 

Pft = -^ + 0(CV), (91) 

which is linear in (tc. This, of course, differs from the Mean Field result. We now argue 
that this naive interpretation is incorrect and that to order (tt we must be more careful in 
defining the reflection and transmission amplitudes. Physically this is due to the fact that 
the angle of the wave has been changed to order ^Tr by the scattering. This change in angle 
should not affect the transmission amplitude. More concretely recall that to the left of the 
scattering region the scattered wave has the form 

v^scat = ^sin(Cvr)x(y)e-^'^^ (92) 

(e is negligible for large k^) and that xiv) is a "sawtooth" with discontinuities whenever 
y=n^. We have discussed previously that x must somehow exponentiate to give the correct 
plane wave in the y direction. It is thus incorrect to read off the transmission coefficient 
naively from Eq. (^2]) which would give 

rft = l + <7rx + 0(CV), (93) 

In fact the contribution of C, is simply to modify the plane wave to take into account the 
"bending" of the particle in the magnetic field and the correct t to this order is one. 
Thus the apparent 0{(7t) correction to the transmission amplitude is best understood as an 
artifact of our approach and the difficulties which the Aharonov-Bohm phases cause when 
one tries to add up the scattered waves. 

Once we realize that r is one to order (tc it is clear that p must vanish to this order. 
This follows from flux conservation. The fact that the pft which we naively calculated is of 
order (it is a direct consequence of the fact that the "naive" rft has a piece which is linear 
in since 

where we have used the fact (from Eq. (|75|) ) that 

x{y) = x{y = -50 - (95) 

It may of course happen that has real pieces quadratic in (n (to which our calculation 
would be insensitive). In this case these should be included in |Tftp. In fact the above 
argument used numerical results and it is thus just valid in an approximate sense. In any 
case it is evident that the apparent discrepancy between the form of p and r in the "flux 
tube" and the Mean Field case is simply an artifact of our approach. 
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IV. DISCUSSION 



In this paper we set out to understand the vahdity of Mean Field Theory when apphed 
to scattering from infinitesimally small flux tubes for which there is no Lorentz force except 
possibly at the location of the flux tubes. We considered the scattering of a plane wave 
representing an incident charged particle from an array of flux tubes with infinite extent 
in the y direction. By computing the current of the wave beyond the scattering region we 
found that the incident beam was bent by the array of flux tubes. When the flux per tube 
is small the angle by which the beam is bent is given precisely by the Mean Field result. 

Caenepeel and MacKenzie have recently shown that no bending occurs when scat- 
tering from a single flux tube unless an additional interaction besides the Aharonov-Bohm 
interaction is present [^. From our calculation one can see that even in the absense of 
an additional interaction there is precisely the correct (MF) amount of bending from an 
array of flux tubes. Although our calculation was done (for essential technical reasons) for 
a regular lattice, it seems very likely on physical grounds that the main results will persist 
for a random array of flux tubes. 

Despite the absense of a Lorentz force for a lattice of flux tubes there is a way to 
understand how such bending can occur classically. Imagine first taking each flux tube 
to be of finite but small size, computing the scattering, and then letting the size become 
very small keeping the mean field, and thus the flux per tube, constant. One finds that the 
charged particle will have to travel an increasingly large distance before scattering but, after 
many interactions, the mean impulse it receives from all the tubes from which it scatters is 
precisely that which it would have received from the mean field. There is however a catch. 
This only works when the flux tubes are larger than some critical size at which the magnetic 
field inside the flux tube becomes so large that the radius of the path of the particle is equal 
to the radius of the flux tube. Beyond this critical size, the average deflection which the 
particle receives becomes too small. It is only Quantum Mechanical scattering to which the 
Mean Field result truly applies when the flux per tube is small. 



V. APPENDIX: COMMENTS ON THE NUMERICAL WORK 

In this appendix we describe several interesting features of our numerical work which 
were not mentioned in the text. 

We begin by discussing the numerical work leading to Eq. (|74D for \&scat- Recall that for 
ky=Q we claimed that our numerical results lead to a scattered wave of the form 

^ e-*'=^(e + 2sin(C7r)x(2/)), (96) 

with 

x{y) ^ -^^y^ + -^kT ^ ^ 

in the region to the left of the flux tubes. The result ( p6D was obtained by summing the 
series in Eq. (p2[) numerically until satisfactory convergence was achieved. For very small 
this involved adding the contributions from up to a million flux tubes in order to obtain 
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a slope in x accurate to about .2%. As one might expect from the fact that the vector 
potential approaches the MF result exponentially (c.f. Eqs. (|33|) - (^) ) the above formula 
becomes valid within a very short distance of the scattering region (i.e. within a distance of 
about 4 * ^). 

For very large compared to C^r, e was negligible. For /c^ closer to ^vr the maximum 
value of |e| was found to be approximately proportional to ^tt and to have an approximate 
dependence. From ( P7D x is also dominated by a dependence for small k^. For all 
cases investigated (we examined C, as low as .0001 and k^ as low as .005) it has been found 
that 

|e| «sin(Cvr)|x|, (98) 

for small k^. 

The growth of x with small has an interesting interpretation. From a MF point of view 
we would expect the particle to be totally reflected for /c^<2|C|7r. From our perturbative 
approach we see that the results become non-perturbative when Ci^x becomes of order 1 i.e. 
when /c^~|C|7r. In a loose sense our approach predicts the correct "critical" point. 

Numerical work was also required for evaluating O defined in Eq. (§3). In order to 
determine O for the case ky=Q the exact forms of fx and fy were used for small kr. For 
larger kr approximate forms (easily derived from (0) and (|58|)) were employed. As in 
the previous calculation it was necessary to sum up to a million tubes. For k^ very large 
compared to 2C,tt and slightly less than 2tt the form quoted for Vty in ( |55| ) was accurate to 
about .2%. Similarly, VL^J^y was "zero" to a magnitude of about .002. 

For the more general case where ky is nonzero, things are much the same. A small 
difference was that some precision was lost in VL^ when ky«kx. The reason is that the 
number which was expected was proportional to ky/ k^ and thus became very small. 

For large momenta it was found that \l'scat and f2 deviated from the forms quoted in 
the main body of the paper due to diffraction effects. Recall that diffraction occurs when 
circular waves "emanating" from neighbouring sites add up coherently. For small momenta 
there is only one angle at which this occurs and only the "forward peak" is observed. For 
higher momenta it becomes possible to have coherence in more than one direction. In effect, 
the transmitted wave becomes the sum of more than one plane wave. As discussed near the 
end of Section III. A, the scattering angle then becomes ill defined. 

One can derive the limit on k^ and ky by using Eq.(^) for the "critical" fluxes. Since 
our numerical approach is only valid near C=0, we shall for convenience set C=0 for this 
derivation. Diffraction then starts when either C^i or is equal to zero so that the limit 
on the momenta is easily seen to be 

{kl + kl)"H+\ky\i<2T:. (99) 
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